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i-^3 ■ The minisuperspace model in 3 + d spatial dimensions with the matter, described by 

the bag model, is considered with the aim to estimate the probability of creation of the 
compactificd extra dimensions in nuclear collisions. The amplitude of transition from 3- 
to 3 + d-dimensional space has been calculated both in the case of completely confined 
^ ■ matter, when contribution of radiation is ignored, and in the case of radiation domi- 

nation, when the bag constant is negligible. It turns out that the number of additional 
' dimensions is limited in the first regime, while this is infinite in the second one. It is 

■ shown that the probability of creation of extra dimensions is finite in the both regimes. 

\J^} • Keywords: Quark-gluon plasma; bag model; Wheeler— DeWitt equation; extra dimen- 

' sions; transition amplitude. 
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1. Introduction 

As it is expected, the Large Hadron Collider (LHC) will be able to detect quantum 
■ gravitational effects. Theoretically, the high-energy particle beam at the LHC can 

cool by boiling off gravitons into the extra space-like dimensions.^ More exotic 
gravitational objects, such as m icro black holes, are assumed to be also produced at 
the TeV energies (see RefsPEEland references therein). In fact the LHC can explore 
these additional dimensions through apparent breaking the energy conservation, as 
well as the appearance and disappearance of gravitons from extra dimensions as 
envisioned by the superstring theory. 

Expectations to observe strongly gravitational effects at the TeV energies mainly 
follow from Ref. 5 , where it is supposed that the electroweak scale ttiew ~ 10 3 GeV 
is the only fundamental short distance scale in nature. If the gravity is unmodified 
down to the Planck length, one might assume that the fundamental mass Mf in 
(3 + g?) + 1-dimensional theory is of ttiew order at the some value of d (the number of 
the compactificd extra dimensions within the framework of Kaluza-Klein scenario). 
As is pointed out in Ref!^, to satisfy this requirement, parameter d should be equal 
to 2. At d = 0, the usual Plank mass, M P i = y/hc/(16irG) « 1.72 • 10 18 GeV, 
is restored. Moreover, such an approach allows one to determine the fundamental 
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scale L of extra dimensions: 



2 



L d ~ _MzL , i , 

Af f 2 + d ' (1) 

Because the inverse square law of gravitational attraction has only been verified 
at distances greater than 1 mm (see RefP and references therein), it is now sug- 
gested that the huge disparity in interaction strengths is caused by the "leaking" of 
gravitational field into the extra dimensions which can be as large as 1 mm with- 
out violating known laboratory, astrophysical, or cosmological data. In this picture, 
the particles and forces are confined to a 3-brane embedded into multi-dimensional 
universe, the only force in the higher-dimensional bulk is gravity. 

Since there is no evidence of the presence of extra dimensions in low-energy 
(classical and quantum) processes, one can assume that an appearance of them is 
induced by the highly dense state of energy and matter like a thermalized quark- 
gluon plasma, predicted by quantum chromodynamics (QCD). Thus, the purpose 
of this paper is to determine the probability of creation of additional dimensions in 
ultrarelativistic heavy- ion collisions. 

The mechanism of creation of extra dimensions can be divided in two steps: 
i) generating the gravitational field in our three spatial dimensions by the matter 
source; ii) transferring the part of energy of this field component into the extra 
dimensions, where new (extra) modes of gravitational field arise. Such a mecha- 
nism is similar to the transfer of the "shaking" (deformation) between slices within 
the elasticity theory. Note that this analogy looks straightforward bec ause t here is 
mathematical relation between general relativity and elasticity theoryf^lSI 

Here we develop the phenomenological model including the matter source, de- 
scribed by the bag model, and the gravitational field, represented by two scale fac- 
tors reflecting two physically different groups of dimensions, the topology of which 
is fixed. The coupling constant (gravitational constant in (3 + d) + 1-dimensional 
space-time) in this model is defined from Eq. (JTJ) . 

Although we will apply the methods and achievements of quantum cosmology, 
but there are several differences between it and our model. Here we assume that 
the 4-dimensional Universe already exists while the quantum cosmology is dealing 
with the creation and evolution of Universe from "nothing" . 



2. Model Formulation 

Let us focus on formulation of the model containing two interacting participants: 
gravitational field and the dense matter being in the quark-gluon phase. As usual, 
we will assume that the homogeneous and isotropic matter lives in 3-dimensional 
space while the gravitational field propagates in 3 + d dimensions. Of course, the 
one time-like evolution parameter is also presented in the model and is common 
for matter and field. We then can imagine that the dense matter, produced in 
nuclear collisions, plays the role of source of gravitational field in three dimensions, 
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whose excitation induces an emergence of the gravitational field modes in extra 
dimensions, the number of which we would like to estimate. 

Let the matter be concentrated in the fireball of radius rf ~ 7 fm, and all of 
three-dimensional phenomena under our consideration are carried out in the bulk 
of this region. The thermodynamic properties (pressure P and energy density e) 
of the quark-gluon phase are described within the framework of the bag modelp^ 
taking into account only light quarks (and gluons), 

P=^~B, + (2) 

where T is the temperature, g — 37 is degeneracy factor, B 1 ' 4 = 200 MeV is the 
bag constant. 

To evaluate the probability of creation of extra spatial dimensions, we should 
obviously to appeal to the quantum picture which provides the appropriate frame- 
work to reach this goal. More precisely, we should deal with the case of quantum 
gravity and/or cosmology in order to find quantum states of (micro) universe and 
to calculate an amplitude of transition between 3- and 3 + (^-dimensional spaces. 

Here we will not account for the topology and geometry variations of extra 
space, which may emerge. Even making use of the gravitational path integral, the 
sum over all topologies is not computable, since there is no algorithm for deciding 
whether two 3 + d-dimcnsional manifolds have the same topology. Moreover, there 
is the problem (even in dimension two) that there exist many more geometries of 
complicated topology than there exist of simple topology. It makes impossible to 
define the theory non-perturbatively in an unambiguous and physically satisfactory 
manner. In higher dimensions these problems are unmanagable. For these reasons 
we are forced to consider the model with the fixed topology of (extra) space, what 
allows us to apply canonical quantization of gravitational field. 

Note that canonical quantization is useful for investigations of the system dy- 
namics, although it requires to impose boundary conditions and to solve the prob- 
lem of the operator ordering. These problems are general for quantum cosmology 
(minisuperspace) models, based on Wheeler-DeWitt equation, and have conceptual 
meaning in the context of the Universe evolution understanding. 

Coming back to the problem of fixing topology, let the extra dimensions be com- 
pactified and described by periodic variables <f> n € (0, 2ttL), where L is the funda- 
mental length of the higher dimensional theory, n = l,d. Then extra <i-dimensional 
manifold £^ is compact and isomorphic to 

SxS...xS . (3) 

d 

Recall that in the original Kaluza-Klein models, space-time was taken to be, at 
least locally, a product manifold of the form R4 x E^, where is some compact 
"internal space" , in the first instance a circle. The effective space-time R4 must thus 
be identified with the quotient of the higher dimensional space-time with respect 
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to S d 

For the sake of simplicity, we can exploit a common scale factor, say 6, for 
additional dimensions. If b = 0, these dimensions are absent, and we concern with 
the theory in four-dimensional space-time. Such an approximation does not permit 
us to investigate the geometry variations for the given topology. Nevertheless, it 
may be effective if the sizes of dimensions are roughly of the same order. 

Similarly, it is possible to introduce a common scale factor for our three dimen- 
sions. In the other words, it means that the 3-space is supposed to be conformally 
flat and isotropic. 

At x 2 + y 2 + z 2 < rf , where the matter is highly dense, the potentials of homo- 
geneous gravitational field, a and 6, are introduced as follows 



d 

ds 2 = -dt 2 + a 2 {t)(dx 2 + dy 2 + dz 2 ) + b 2 (t)^d(j) 2 nl (4) 



where we use the gauge-fixing condition (N — 1, where N is the lapse function) 
resulting in the formalism with "Hubble time" . 

On the other hand, the metric becomes asymptotically Minkowskian and d = 
(and b = 0) at x 2 + y 2 + z 2 > r 2 . 

The quantum states of this model are described by wave function (a, 6|3 + d) 
dependent on d. Therefore the amplitude of transition "3" — > "3 + d" at the fixed 
6, playing role of parameter, is defined as 



dfi(a)(3\a,0) (a,b\3 + d) 



(5) 



where /i(a) is the integration measure. 

Having got the normalized wave function (a, 6|3 + d), we can also find more 
probable values of b at the fixed a and d from the probability density profile, 

\(a,b\3 + d)\ 2 . (6) 

Another possibility to describe extra dimensions consists in the finding of the 
wave function (3 + d\a, 61,62, ...,bd) dependent on 1 + d scale factors. Such an 
approach could be used to study the different configurations (geometries) of 3 + d- 
dimensional space. However, we limit ourselves here by the simple case depicted 
before. 



3. Hamiltonian Dynamics 

In this Section, starting from an action of relativistic system "gravitational 
field+ideal fluid" , we reformulate the dynamics in the terms of canonical variables 
needed for further canonical quantization. 

a Sometimes, authors identify our space-time with a submanifold rather than a quotient. 
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As mentioned in previous Section, we assume that the matter is located in the 
region of 3-dimensional space at x 2 + y 2 + z 2 < r 2 . Therefore there is the following 
energy conservation in three spatial dimensions: 

i + 3-(e + P) = 0, (7) 
a 

where the dot means derivative with respect to t. 

Since the thermodynamic properties of the matter are described by Eqs. ([2|, 
resulting in equation of state, 

p=r-- 3 B, ( 8) 

integrating Eq. ([7]) , one obtains that the dependence of the energy density e on the 
scale factor is 

e = B+^, (9) 
a 1 

where £ — 7r 2 gT 4 /30, and the scale factor at t = to is assumed to be do = 1. 

Thus, the action functional of the (quark-gluon) ideal fluid in the rest frame of 
the matter is of the form: 

Su = — I ea i dtd i x 



= -Vf J (Ba 3 + ^ dt, (10) 

where V{ = Airrf /3 is the fireball volume. 

Note that the successive d eriva tion of Sm can be done, for instance, by means 
of formalism elaborated in RefP] 

To write down the action functional for the gravitational field, we have to find 
the scalar curvature of (3 + d) + 1-dimensional space-time. The term containing the 
four-dimensional curvature is often missing from the action. However, in general, 
this term seems to be essential since it is generated as a quantum correction to 
the matter action. Note that this quantum correction typically involves an infinite 
number of terms of higher order in curvature. Nevertheless, these terms are com- 
pletely suppressed in our consideration and the action functional is formulated in 
accordance with Einstein. 

The non- vanishing Cristoffel symbols of second kind are 

r° — r° - r° - ™ r° - hh 

1 xx — 1 yy — 1 22 — 1 0„0„ — 00 ' 

v x — v v — v z — - r^" — - mi 

1 Qx — 1 Qy — 1 02 — 1 i 0*„ _ J' y 11 ) 

and the Ricci scalar is presented in the following form: 

2 d 2 (n 3 k d ) 

EM] = GM] + _^_J, (12 ) 
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where one can ignore the last term in r.h.s., resulting in total derivative with respect 
to evolution parameter. We have introduced the notation: 



'- 2 Mb „, ^b 2 



G[a,b] = -6^-6d-^-d{d-l)^. (13) 
ab tr 



a 



As was discussed above, the gravitational field is supposed to be homogeneous 
and non-trivial inside three-dimensional fireball while the metric in the space be- 
yond fireball tends asymptotically to the usual Minkowskian one (without additional 
dimensions). Therefore it is possible to limit the region of integration in the action 
functional by the volume filled by the quark-gluon matter. Such an approximation 
means that we are neglecting by the surface contribution, where the gradient of 
pressure, gravitational potentials and curvature, in principal, can be large enough. 
This approximation is valid, if the profiles of matter distribution and scale factors 
are preserved. It is easy to understand that the matter distribution changes are 
related with longitudinal hydrodynamic waves, propagating with the sound speed 
c s = J dPj de s=s 0.57c. Then the characteristic time of these changes can be eval- 
uated as r c = rf/c s w 12.3 fm/c and is approximately equal to the life-time of 
fireball. Therefore we conclude the time of process must be limited, t c 3> At. 

Thus the Lagrangian function of the system, which we consider, is written as 

L = V f [M^G[a, b]a 3 b d - U(a)] , (14) 

where we have applied the relation M 2+d = Mpj J d d (f> n and introduced the matter 
potential U(a) = Ba 3 +£/a. 

Canonical momenta conjugated to scale factors a and b respectively are 

Pa = ^ = -6M^V { a 2 b d [ 2- + dl ) , (15) 

Pb = ^- = -2dM| 1 \4a 3 6 d - 1 ( (d - l)^ + 3^ ) . (16) 

Resolving these relations with respect to velocities d, b, one finds that 

a = 2(2 + d)M^V f ab d {— Pa a Pb ) ' {U) 
1 fb 




Canonical Hamiltonian of the system is 

+ VfU(a). (19) 



H = dp a + bp b - L 

_ d(d - l)a 2 p 2 a - 6dabp aPb + &b 2 p 2 b 



12d(d + 2)M 2 l V { a 3 b d 

Introducing the operators instead of classical variables, we will use this function for 
formulating Whcclcr-DcWitt equation in the next Section. 
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4. Quantization. Probability of Creation of Extra Dimensions 

Before quantization is performed, let us make some modifications. 

Due to general covariance, the system at the classical level should be described 
by Hamiltonian constraint H — 0, which can be re- written as 

2 



12 



bpb 



l -ap a } = a 2 [p 2 a - >c 2 ab d V(a)\ , 



(20) 



d{d + 2) 2 

where we have introduced the coefficient x 2 = 24M^V 2 / (he) 6 w 2.5 • 10 36 MeV~ 4 . 

Since there is no correct limit for Hamiltonian constraint at d — > 0, we perform 
a canonical transformation. It is convenient to introduce new variables £<j = ab d ^ 2 , 
Vd = b d (dependent on d) instead of scale factors a, b. The conjugate momenta are 
related as 



Pi 



= Pab- d / 2 , 



Pa 



VdPi, 



p b = dv\ 



-1/d 



2v d 



Pi 



In the terms of new canonical variables, we obtain 



12d 
d + 2 



2 2 
V dPv 



2 2 v 



-JJ 



At d = 0, this equation is reduced to the following one: 



2 2 Co jj 

Pt-H —=U 



= o, 



(21) 
(22) 

(23) 

(24) 



where £o = a , v a — 1> Pi = Pa as must be. 

Making use of the scheme of canonical quantization, one has to replace momen- 
tum variables with differential operators: 



Pi 



d d 

-«tt^, Pv — ► —l- 



(25) 



d£d dv d 

The Wheeler-DeWitt equation with the normal ordering of operators (which is 
usual for quantum field theory) is written as follows 



I2d o d 2 



d + 2 



dv 2 



5d uu d y/V d \^/V di 

At d = 0, it is reduced to the ordinary differential equation: 



*o = 0. 



(26) 



(27) 



Note that this connection makes sense due to the form of (space-time interval 
and) Hamiltonian where parameter d can be of fractal (non-integer) dimension. 
Moreover, if it is possible to change d continuously, the topology ^ of extra di- 
mensions becomes also fractal. However, we will not discuss here aspects of the 
fractal topology construction. 
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In some cases of U{x), these equations can be solved exactly. However, since 
k 2 U(x) is quite large in our consideration, the quasiclassical WKB approximation 
is applicable enough. 

Note also that, we are not imposing boundary conditions in our minisuperspace 
model. Hence such a problem formulation can be treated as pseudo-Euclidean anal- 
ogy of the Hartle-Hawking "no boundary" proposal. 

In the next subsections we will investigate the particular regimes of the bag 
model equation of state. 



4.1. The Case of Confined Matter 

Let B ^> £ and, therefore, U(a) « Ba 3 . In this case (of cold matter), the matter 
pressure is negative and is approximately equal to — B. The dynamics of gravita- 
tional field is strongly damped in this regime. Just decreasing B, a sufficient role of 
the field can be expected. However, such a (confined) state of matter, predicted by 
QCD, should be investigated for better understanding of phenomena, which emerge 
under unification of strong and gravitational interactions. 

Now let us describe briefly the general scheme used in further calculations. 
Limiting ourselves by a first quantum correction, the wave functions are found in 
the form: 



# ~ e lS/h , S^S + -S, 



qi 



(28) 



where we write down explicitly Plank constant h; So is the classical action, which 
satisfies the Hamilton- Jacobi equation: 



12d 
d + 2' 



dSo 
dv d 



e d 



dSo 



II 



The first quantum correction S q is found from equation: 

dS dS q / . 2 d 2 S a I2d f^ 2 dS dS q 2 d 2 S 

+ v 



d + 2 



dv d dv d 



Substituting U(a), the wave functions in WKB approximation are 

v d 



■ exp in 



B2 



3 6 — d Vd 



(29) 



(30) 



(31) 



B 



■ exp ix- 



-a" = 



exp iK- 



3/2 



(32) 



Looking at ^d, we see that the number of extra dimensions is limited by d — 
6. At d > 6 the wave function becomes real and exponentially decreasing. Since 
K\f~B k, 6.3 • 10 23 , ^d>6 vanishes. On the other hand, the wave function 



Vd 

— exp 

t,d 



-IK 



B 2 + d Q 
3 6 — d Vd 



(33) 
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becomes exponentially increasing at d > 6 and, therefore, this state is unphysical. 
It confirms that the number of extra dimensions in this model does not exceed six. 

Using the wave functions found, let us consider the transition amplitude from 
3- to 3 + e?-dimensional space at the fixed Vd- We assume that the all compact 
dimensions are opened i.e., v c i = 1. In this case, the wave functions are described 
by a common expression: 



&|3 + i) = ^exp^|^j. 
This wave function is normalized with respect to the integration measure: 



\(U^ + d}\ 2 dU = h 



(34) 



(35) 



where the lower limit is chosen to be H > 0. If S = 1, one has exactly the scale factor 
of Minkowski space-time. It means that the integration over the region ^ € (EE, oo) 
accounts for the expanding geometry only. The process of collapsing (when ^ G 
(0, EE)) should be ignored in the system with negative pressure in order to eliminate 
unphysical situation. Mathematically, the integration in the range £d € (0, EE) brings 
a divergent term. 

The transition amplitude is defined as 



T(d) 



(3|&)(&|3 + d)d& 



The result of computations is presented as 



I +%zU[ 1,-, 



T(d) = 

where U(a, b, z) is the Kummer function 




(36) 



(37) 



(38) 



At z — > oo and d > 0, the asymptotic expression for the transition amplitude is 
T w l/(3z). At d -> (z -> 0), we have to use Eq. (|37|) and obtain T(0) = 1. 

Since the number of extra dimensions is limited, we can construct the continuous 
probability density as 



V(d) = T 2 {d) / J T 2 { 



x)dx 



(39) 



which depends on the model parameter fi = ><\ABEE 3 /3. Here we will not explore 
the question: what happens if we put d < 0. 

The behavior of V(d) is demonstrated in Fig. 1. We may immediately conclude 
that the probability of creation of few extra dimensions is finite] 
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Fig. 1. Probability density of creation of extra dimensions at fi = 0.11 (solid line), 0.01 (dashed 



Now, making use of the probability density V(d), it is interesting to investigate 
the average number of additional dimensions: 

6 



This is function of parameter fi. Thus, we get that (d) — > 3, when /i — * 0, and 
(d) — ► at /i — ► oo. In order to reproduce (d) — 2, what is predicted by ADD 
paradigm, one has /i « 0.11. In principal, we should treat A c ff = BE 6 as the 
vacuum energy density of QCD in the present regime. To realize (d) = 2, we get 
that A e ff « 4.36 • 10 -38 MeV 4 . Since the cosmological constant is often interpreted 
as the vacuum energy , one can conclude that this quantity A = 3.956 • 10 -35 MeV 4 



4.2. The Case of Radiation Domination 

Now let us assume £ 3> B and U(a) ps £/a. It means that the temperature of the 
quark-gluon gas is much higher than B 1 / 4 = 200 MeV. This case is an interesting 
one because the energy density of matter (and the matter temperature) can be 
changed due to gravitational effects like an origin of additional dimensions, while 
the bag constant is supposed to be unchanged. 

Since the energy density £ is the function of the temperature T, which may be 
changed, let us introduce quasi-momentum k = >c\J~£~ as the free parameter and the 
characteristic quasi-momentum K ps 2.2 • 10 23 for T = 200 MeV. In principal, k 
should be bigger than K in the regime under consideration. 

The Wheeler-DeWitt equations without boundary conditions for the system, 
where radiation dominates, are solved exactly and give us the wave functions: 



line). 




(40) 



o 





(41) 
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, '- 4 MWm w,) ' (42) 

where is the Hankel function; A is arbitrary number such that 3?A > 0. 

Note that is infinite at £d — 0, if A(A — 1) > 0. Indeed, the further analysis 
is proving that A(A — 1) < and A = 1/2 — ifj,, /j, S M. In the last case, the index 
of the Hankel function is purely imaginary number, resulting in vanishing the wave 
function at £<j = 0. 

Remembering that k > 10 23 , let us consider an asymptotic expression (up to 
constant phase) for ^d- 

~ ex P (l + *A(A -D^4+0 ( (^y) ) - (43) 

where we will omit the term of order O (l/(fc 2 ^)). 

On the other hand, the WKB approach leads to the following classical action 
and the first quantum correction: 



go = ^ + ^fc 2 -v / ^ln(2 C " 1 v ^ ' Sd ' v I -^^linw (44) 



where c, Ci, C2 are arbitrary constants. 

Using these formulas, we try to analyze the consistency of the approximate 
solutions and to connect the constants. 

First of all, let us expand Sq and S q in the same manner like vp^: 



«, = ft - ^In (2WH) - f v/^ta* + - ^ + O (^) , (46) 

S «- C =4"-^ 1 "< 2 ^)-^i + °(4)' <47) 
Comparing exp (iSo + S q ) with from Eq. (|43|) . it is convenient to take 

d = ^ln(2fcv^), C 2 = ilnfc + y^ln(2^), c=^ S 2 , (48) 

where s is regarded as a new parameter instead of c. 

Comparing this expression with Eq. (|43|) again, we derive that s = 2zA and 
A(A — 1) = — (s 2 — 2is)/4. In order to be s 2 — 2is € R (the Hamiltonian operator 
has to be (quasi)Hcrmitian), one finds parametrization, s = 2/i + i, resulting in 
s 2 - 2is = 1 + 4/i 2 , where 

We interpret the case of ft, = as the ground state of the system. Without 
loosing generality, we will focus on this situation. It is clear that the excited states 
are less probable than the case at a hand. 
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Now let us discuss the normalization condition for the wave functions. We will 
assume again that vj, = 1. Then the scale factor £<j can run from to oo in this case 
because there is no problem with the negative pressure and collapse. Because of 
I2 symmetry of the space-time interval (with respect to the sign inversion of scale 
factors), the range of £<j can be further extended to (—00, +00). 

Note that the asymptotic expression for ^d has the pole at £d = 0. This property 
is important for normalization procedure. To determine the integration contour near 

= 0, we appeal to the exact solution (|42|) : the symmetry of \^>d\ 2 profile in the 
ranges £<j S (—00, 0) and £d € (0, +00) is preserved only if £<j — > £<j + i0 + . 

Using well-known formula for Dirac <5-function: 



e lkx dx = 5(k), (49) 



1 

2^ 

the wave function |3, k) = exp (ik^d) / \/27r at d — is normalized as 

(3,q\3,k) = 5(k-q). (50) 

This relation reflects simply the energy conservation in three dimensions. 
The normalization condition for the wave function at d > 0, 

|3 + d ' fc) = ^ eXP( ^ ) ( 1 - i 2(l^)4)' (51) 



looks like 



00 

(3 + d,q\3 + d,k) w -L /" d^e 



1 



3d /l 1\ 1 



2(rf + 2) 9/ +»0 



(52) 



We see that the energy density of matter can be varied in the presence of extra 
dimensions! 

Note that, taking term 0(l/(fc 2 ^)) into account in |3 + d, k), one gets term of 
0((k — q) 2 / (k 2 q 2 )) order after integration. Therefore, our approximation is valid for 
the small energy deviation \k — q\/k. The 0-function is defined such that 8(0) = 1/2. 

Having got all necessary ingredients, we are able to calculate the amplitude of 
transition from 3- to 3 + d-dimensional space. This quantity is 

(3, g|3 + d, k) = 5(k -q)- ^—±6(q k). (54) 

The Heaviside ^-function in the last term corresponds to the cooling of the matter. 
It is understandable that this process is related with the transferring of energy from 
the matter to the extra dimensions. 

Moreover, there is another effect consisting in the changes of the system volume. 
If Vd 7^ 1, we have to replace q with q/ ^/vd in the last formula. 



I 



I 



June 22, 2008 16:55 WSPC/INSTRUCTION FILE Probability 



Probability of creation of extra dimensions... 13 

The probability of creation of extra dimensions in the case of radiation domi- 
nation is defined as 

oo 

W = J \(3,q\3 + d,k)\ 2 dk. (55) 



As usualpl W e should make the following replacement: 

L 

5 2 {k ^ q )^5{k-q)^ L J e^-^dtd = 5(k - q). (56) 
-L 

Limiting by the term of order 1/q, one finds 

W = l-^T i - (57) 
2(d + 2)q V ' 

This formula explains why the high energy (or temperature T) is needed to create 
extra dimensions. However, the model predicts the infinite number of additional 
dimensions in the contrast with the case of the confined matter. Since 

W^oo = 1 - (58) 
2q 

it is impossible to find the probability density in our approximation. 



5. Conclusions 

In this paper the minisuperspace model in 3 + d spatial dimensions with the matter, 
described by the bag model, has been considered with the aim to estimate the 
probability of the origin and the number of the compactified extra dimensions in 
nuclear collisions. The gravitational constant in (3 + d) + 1 space-time is introduced 
on the base of ADD paradigm^ and the topology of extra dimensions is supposed 
to be toroidal in accordance with Kaluza-Klein scenario. 

For the sake of simplicity, the amplitude of transition from 3- to 3+c?-dimensional 
space has been calculated in two regimes of the bag model: i) in the case of com- 
pletely confined matter (what is equal to the model with the cosmological constant); 
ii) in the case of radiation domination, when the bag constant is negligible. 

In the first case, the number of additional dimensions is limited by d = 6, that 
is agreed with superstring theory prediction. 

It is shown in the second case that the energy density (the temperature) of the 
matter is changed in the presence of extra dimensions, the number of which is not 
limited in approximation used. 

It turns out that it is possible to evaluate the mean value of extra dimensions 
in the first case only and to obtain (d) = 2 in agreement with ADD predictions. 

It is important to emphasize that the probability of creation of extra dimen- 
sions is finite in the both regimes. The probability for intermediate regimes is not 
computed here and should be investigated. 
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